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Three techniques for performing gauge-invariant, noncompact lattice simulations of nonabelian gauge theories 
are discussed. 

In the first method, the action is not itself gauge invariant, but a kind of lattice gauge invariance is restored 
by random compact gauge transformations during the successive sweeps of the simulation. This method has been 
applied to pure 5(7(2) gauge theory on a 12'* lattice, and Wilson loops have been measured at strong coupling, 
/3 = 0.5. These Wilson loops display a confinement signal not seen in simulations performed earlier with the same 
action but without the random gauge transformations. 

In the second method, the action is gauge symmetrized by integrations over the group manifold. 

The third method is based upon a new, noncompact form of the action that is exactly invariant under lattice 
gauge transformations. The action is a natural discretization of the classical Yang-Mills action. 



1. INTRODUCTION 

In Wilson's formulation [Q of lattice gauge the- 
ory, the basic variables are the elements of the 
gauge group, not the fields of the theory. Using 
this formulation Creutz |^ in 1980 displayed ev- 
idence for confinement in both abelian and non- 
abelian gauge theories. 

The Wilson action is exactly invariant under 
lattice gauge transformations and reduces to the 
classical action when the group elements are close 
to the identity. Wilson's elegant formulation is 
ideal at weak coupling. 

But at strong coupling, the group elements are 
often far from the identity, and the Wilson action 
is quite different from the classical Yang-Mills ac- 
tion. It is unknown whether Wilson's formulation 
is accurate at the moderately strong couplings 
where it is really needed. Among the symptoms 
of trouble are the false vacua |^ of the Wilson 
action which at strong coupling affect the string 
tension and the confinement signals exhib- 
ited by Creutz |Q in his simulations of abelian 
gauge theories. 

To examine these questions, some physicists 
have introduced lattice actions that are noncom- 
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pact discretizations of the continuum action with 
the fields as the basic variables. For U{1) these 
noncompact formulations are accurate for all cou- 
phng strengths |^ ; for SU (2) they agree well with 
perturbation theory at very weak coupling |Q . 

Patrascioiu, Seller, Stamatescu, Wolff, and 
Zwanziger performed the first noncompact sim- 
ulations of SU{2) by discrctizing the classical 
action and fixing the gauge They saw a 

Coulomb force. 

Later simulations |^ were carried out with an 
action free of spurious zero modes, for which it 
was not necessary to fix the gauge. The Wilson 
loops of these simulations showed no sign of quark 
confinement. A likely explanation of this negative 
result is that these noncompact actions lack an 
exact lattice gauge invariance. 

The first gauge-invariant noncompact simula- 
tions were carried out by Palumbo, Polikarpov, 
and Veselov [|l^ and were based on earlier work 
by Palumbo et al. ijlTI . They saw a confinement 
signal |ic|] . Their action contains five terms, con- 
structed from two invariants, and involves auxil- 
iary fields and an adjustable parameter. 

The present paper discusses three ways of per- 
forming gauge-invariant noncompact simulations 
with an action that is a discretization of the clas- 
sical Yang-Mills action without extra terms or 
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fields. 

In the first technique, one subjects the fields 
to random compact gauge transformations during 
the Monte-Carlo updates The random gauge 
transformations restore a semblance of gauge in- 
variance |l2| . This method has been used to mea- 
sure Wilson loops at strong couphng, (j = A/ g'^ = 
0.5, on a 12^ lattice in a noncompact simulation 
of SU{2) gauge theory without gauge fixing or 
fermions. In this simulation the Wilson loops fall 
off exponentially with the area of the loop at a 
rate that is over two orders of magnitude greater 
than the negligible rate seen in an earlier simula- 
tion with the same action but without the random 
gauge transformations. 

The second technique is to symmetrize the ac- 
tion itself by subjecting it to a full set of gauge 
transformations. Here one inserts a gauge trans- 
formation at each vertex and integrates over the 
gauge group. This technique is simpler to perform 
if the gauge transformations are noncompact. 

The third method is based upon a new noncom- 
pact action that is exactly invariant under lattice 
gauge transformations. The action is a natural 
discretization of the classical Yang-Mills action. 
Three ways of interpreting this action are out- 
lined. 

2. THE FIRST METHOD 

2.1. The Action 

In both the earlier simulations without gauge 
invariance and the present simulation with ran- 
dom gauge transformations, the fields are con- 
stant on the links of length a, the lattice spacing, 
but are interpolated linearly throughout the pla- 
quettes. In the plaquette with vertices n, n -|- e^, 



Cj/ , and n 



e^, the field is 



+ K + l-^)A;(n) 



(1) 



and the field strength is 

+ gf^,Al{x)Al{x). (2) 

The action S is the sum over all plaquettes of the 
integral over each plaquette of the squared field 



strength, 

^ = Ey j dx.dx^F^^Mf 



(3) 



The mean value in the vacuum of a euclidean- 
time-ordered operator Q is approximated by a 
ratio of multiple integrals over the ^°(n)'s 

which one may compute numerically. Macsyma 
was used to write most of the Fortran code ||l^ 
for the present simulation. 

2.2. Gauge Invariance 

To restore gauge invariance, the fields are sub- 
jected to random compact gauge transformations 
during every sweep, except those devoted exclu- 
sively to measurements. At each vertex n a ran- 
dom number r is generated uniformly on the in- 
terval (0, 1); and if r is less than a fixed proba- 
bility, set equal to 0.5 in this work, then a ran- 
dom group element U (n) is picked from the group 
SU{2). The fields on the four finks coming out of 
the vertex n are then subjected to the compact 
gauge transformation 

^^igaA'l^{n)Ti ^ 6^*^"'^^ [/(n)''' (5) 

and those on the links entering the vertex to the 
transformation 

g-igayl^''(n-e^)T6 ^ [7(^)6"'^"^°' . (6) 

At weak coupling one should use random gauge 
transformations that are suitably close to the 
identity. 

2.3. Wilson Loops 

The quantity normally used to study confine- 
ment in quarkless gauge theories is the Wilson 
loop W{r, t) which is the mean value in the vac- 
uum of the trace of a path-and-time-ordered ex- 
ponential of a line integral of the connection 
around an r x f rectangle 

W{r,t) = (1/d) (tr7'Te-*»i'^°'^"'^"")o (7) 

where d is the dimension of the generators Ta- 
Although Wilson loops vanish in the exact the- 
ory n4], Creutz ratios x(r, t) of Wilson loops de- 
fined [|15| as double differences of logarithms of 
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Wilson loops are finite. For large t, x(r, t) ap- 
proximates (a^ times) the force between a quark 
and an antiquark separated by the distance r. 

In this simulation the data are not yet suffi- 
cient to allow one to determine the Creutz ratios 
beyond the 3x4 loop. The Wilson loops there- 
fore have been fitted to an expression involving 
Coulomb, perimeter, scale, and area terms. 

2.4. Measurements and Results 

It will be useful to compare this simulation with 
an earlier one ||^ in which the fields were not sub- 
jected to random gauge transformations. Both 
simulations were done on a 12'* periodic lattice 
with a heat bath. The earlier simulation consisted 
of 20 independent runs with cold starts. The first 
run had 25,000 thermalizing sweeps at inverse 
couphng P = 2 followed by 5000 at /? = 0.5; the 
other nineteen runs began at /9 = 0.5 with 20,000 
thermalizing sweeps. There were 59,640 Parisi- 
assisted jl^ measurements, 20 sweeps apart. 

The present simulation with random gauge 
transformations is very noisy. It consists of 21 
runs of which 20 were from cold starts and one 
from a random start. Each began with 20,000 
thermalizing sweeps. After the first few hun- 
dred thermalizing sweeps in each run, the average 
value of the action was stable. 

Wilson loops have been measured every five 
sweeps for a total of 1,313,202 measurements. 
The values of the Wilson loops so obtained are 
listed in the table. The errors have been esti- 
mated by the jackknife method, with all measure- 
ments in bins of 4000 considered to be indepen- 
dent. 

The Wilson loops of the gauge-invariant sim- 
ulation fall off much faster with increasing loop 
size than do those of the earlier simulation. Be- 
cause the data do not accurately determine all 
the Creutz ratios, I have fitted both sets of loops, 
including the non-diagonal loops, to the formula 

W{r,t) « Q'i+bit/r+r/t)-2c(r+t)-drt j^g^ 

in which a is a scale factor, b a Coulomb term, c a 
perimeter term, and d an area term. For the sim- 
ulation without random gauge transformations, I 
found a « 0.26, b w 0.20, c « 0.39, and d w 0.00. 
For the simulation with random gauge transfor- 
mations, I found a « 0.57, b w 0.15, c w 0.51, and 



d w 0.18. In the gauge-invariant simulation, the 
coefficient of the area-law term is over two orders 
of magnitude larger than in the earlier simulation 
which lacked gauge invariance. 



Noncompact Wilson loops at (3 = 0.5 
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Invariant 


1 X 1 


0.402330(6) 


0.254558(5) 


1x2 


0.208096(5) 


0.082528(6) 


2x2 


0.085426(4) 


0.018709(4) 


1x3 


0.111792(5) 


0.027715(5) 


2x3 


0.040008(3) 


0.005224(3) 


3x3 


0.018080(2) 


0.001431(5) 


1x4 


0.060451(5) 


0.009352(3) 


2x4 


0.019212(2) 


0.001502(5) 


3x4 


0.008517(1) 


0.000407(3) 


4x4 


0.003993(1) 


0.000121(5) 


1x5 


0.032726(5) 


0.003152(3) 


2x5 


0.009269(1) 


0.000434(5) 


3x5 


0.004041(1) 


0.000118(3) 


4x5 


0.001889(1) 


0.000039(5) 


5x5 


0.000893(1) 


0.000013(3) 


1x6 


0.017717(4) 


0.001064(5) 


2x6 


0.004474(1) 


0.000123(1) 


3x6 


0.001919(1) 


0.000033(1) 


4x6 


0.000896(1) 


0.000009(1) 


5x6 


0.000423(0) 


0.000003(1) 


6x6 


0.000201(0) 


0.000003(1) 



To exhibit the renormalized quark- antiquark 
potential, I have plotted in the figure the negative 
logarithms - logio \^e-'^(r/t+t/r)+2c(r+t)y^^^^^y^ 

the Wilson loops with the Coulomb and perimeter 
terms removed. Apart from the uncertain value 
of W{&, 6), the loops of the gauge-invariant simu- 
lation, represented by bullets, seem to display an 
area law; whereas the larger loops of the earlier 
simulation, represented by circles, show an essen- 
tially flat potential. 




1 1 1 1 1 1 r 

5 10 15 20 25 30 35 40 

Area 



Figure 1. The negative logarithms of Wilson 
loops with the Coulomb and perimeter factors 
canceled are plotted against the area rt of the 
loop. The loops of the gauge-invariant simula- 
tion are represented by bullets; those of the ear- 
lier simulation by circles. 



3. THE SECOND METHOD 

While the random gauge transformations of the 
first method can restore a measure of gauge in- 
variance, it clearly would be better to use an ac- 
tion that is itself exactly gauge invariant or has 
been made so. In the second method of perform- 
ing gauge-invariant noncompact simulations, one 
symmetrizes an action that is not itself gauge in- 
variant. To make the action gauge invariant, one 
inserts a general gauge transformation at each 
vertex of the lattice and then independently in- 
tegrates over the group manifold at each vertex 
using the Haar measure. 

One may take any suitable noncompact action 
as a starting point, for instance the action used 
in the first method. Since the action is positive 
definite, the result of the integration over the pa- 
rameters of the group will not be zero. In practice 
one need only integrate over the vertices of the 
plaquettes that contain the arbitrary link that is 
to be updated, fourteen in the case of the action 
(^). The relation between the gauge field before 



and after the gauge transformation is transcen- 
dental in the case of compact gauge transforma- 
tions and affine in the case of noncompact gauge 
transformations. Obviously the latter kind would 
be easier to implement. 

4. THE THIRD METHOD 

The third method of performing gauge- 
invariant noncompact simulations is to start with 
a noncompact lattice action that is intrinsically 
gauge invariant. 

4.1. Gauge Invariance 

Let us first decide what constitutes a gauge 
transformation in this third method. A good way 
to do that is to start with the fermionic action 
density which in the continuum is 



(9) 



A suitable discretization of the free part of that 
action density is 



(10) 



in which n is a four- vector of integers representing 
an arbitrary vertex of the lattice, is a unit 
vector in the /zth direction, and a is the lattice 
spacing. The product of Fermi fields at the same 
point is gauge invariant as it stands. The other 
product of Fermi fields becomes gauge invariant 
if we insert a matrix Ap^in) of gauge fields 

^Vi(n)7^[(l -I- igaAf,{n))tP{n + e^) - (11) 

that transforms under a gauge transformation 
represented by the group elements U{n) and 
U{n + e^) in such a way that 

l+iagA'^{n) = [/(n)[H-m.g^^(n)]f/-i(n+e^).(12) 

The required behavior is 

A'^in) = U{n)Af,in)U-\n + e^) 

+ —Uin)[U-\n)^U-\n + e^)]. (13) 
ag 

Under such a gauge transformation, the usual 
gauge-field matrix Afj_{n) — TaA1^{n) remains in 
the Lie algebra only to first order in the lattice 
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spacing a. Three ways of coping with this prob- 
lem are outlined in Sec. 4.2. 

Let us define the lattice field strength F^^{n) 

as 



+ i9[K{n)A^{n + e^) 
- A^,{n)A^{n + e^,)] 



(14) 



which reduces to the continuum Yang-Mills field 
strength in the limit a — > 0. Under the afore- 
mentioned gauge transformation (p^, this field 
strength transforms as 

K'^iri) = U{n)F^,{n)U-\n + e^ + e,). (15) 

The field strength F^^, {n) is antisymmetric in the 
indices /x and i/, but it is not hermitian. To make 
a positive plaquette action density, we write 



^Tr[Ft,(n)F^,(n)], 



4fc 



(16) 



in which it is assumed that the generators Tq of 
the gauge group have been orthonormalized as 
Tr(raTb) = kSab- Because F^j/(n) transforms co- 
variantly (|l^), this action density is exactly in- 
variant under the noncompact gauge transforma- 
tion (|f|). 

4.2. Three Interpretations 

In general the gauge transformation (^3|) maps 
the usual matrix of gauge fields A^{n) = TaA^(n) 
into a matrix that lies outside the Lie algebra of 
the gauge group. I shall now outline three re- 
sponses to this problem. 

The first response is to note that for group el- 
ements U (n) of the form 



U{n) ^ 6"'"^'"°^" 



(17) 



the gauge transformation (13) to lowest (zeroth) 
order in the lattice spacing a, does keep the 
gauge-field matrix in the Lie algebra. Thus one 
may use the usual matrix A^{n) — TaA'^^{n) of 
gauge fields and retain invariance under infinites- 
imal gauge transformations. 

The second response is to accept the fact that 
the gauge-field matrix A^{n) will be mapped by 



the gauge transformation (|l^) into a matrix that 
lies outside the Lie algebra of the gauge group and 
to use this more-general matrix in the simulation. 
Thus one may use the action (|l^) in which the 
field strength (|lj) is defined in terms of gauge- 
field matrices that are of the more general form 



W-^) (18) 



A^{n) = VAl{n)W-^ + —V {V- 

ag 

where A'^(n) is a matrix of gauge fields defined 
in the usual way, A^^{n) = TaA'^'^{n). Here the 
group elements V and W associated with the 
gauge field Afj,{n) are unrelated to those associ- 
ated with the neighboring gauge fields Ai^{n+e^), 
Ai,(n), and A^{n + e^). I intend to test this 
method in the near future. 

The third response is to take a cue from the 
transformation rule (|l^) and to represent the 
quantity 1 -|- igaA^in) as a element L^{n) of the 
gauge group. In this case the matrix A^(n) of 
gauge fields is related to the link L^in) by 



iga 



(19) 



and the action (16) defined in terms of the field 
strength (^4|) with gauge-field matrix ( |l9|) is, 
mirabile dictu, Wilson's action 



S 



k - ?SiTiLf,{n)Liy{n + e^)Lj,(n + ey)Ll{n) 
2a^g'^k 



(20) 
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